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We study the thermodynamics of the electroweak theory at a finite lepton number density. The 
phase diagram of the theory is calculated by relating the full 4-dimensional theory to a 3-dimensional 
effective theory which has been previously solved using nonperturbative methods. It is seen that 
the critical temperature increases and the value of the Higgs boson mass at which the first order 
phase transition line ends decreases with increasing leptonic chemical potential. 
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I. INTRODUCTION 

The complete thermodynamic description of the elec- 
troweak theory 1 depends on only five intensive variables; 
the temperature T of the system, the strength of the ex- 
ternal U(l) magnetic field Hy and the leptonic chemical 
potentials /ii^. The most studied case is that when only 
the temperature T and the conjugated variable entropy 
are nonzero as it was early understood that at high tem- 
peratures the symmetry of the electroweak theory would 
be restored Q. Much work was devoted to thisproblem 
using, for example, perturbative 1-loop 0, 0, 3 and 2- 
loop [5j effective potential calculations. These only work 
for small scalar self-couplings or for small Higgs masses 
and the full solution of the problem required first a per- 
turbative matching of the full 4-dimcnsional theory to 
an effective 3-dimensional theory [|| . The phase diagram 
of the effective theory was then numerically solved with 
lattice Monte Carlo techniques Q (the phase diagram 
has been studied numerically also with the 4-dimensional 
theory, see [jj). The result is that the phase diagram 
contains a first-order line which ends in a 2nd order crit- 
ical point of Ising universality class Q. Similar tech- 
niques were then applied to solve the phase diagram when 
also Hy and the conjugate extensive variable VBy were 
nonzero 0] ■ The purpose of this paper is to study the 
remaining case; how the phase diagram depends on finite 
chemical potentials related to lepton and baryon num- 
bers and on the conjugate extensive variables, net lepton 
and baryon number densities. 

Thermodynamical properties of the electroweak theory 
at nonzero lepton number density are interesting from 
many points of view. Theoretically, the minimal stan- 
dard model describes Nature to very high accuracy and 
thus it is important that we know the theory completely. 
Especially the partition function is a fundamental con- 
cept and to know it under most general circumstances 
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1 To fix the terminology, in this paper electroweak theory means 
the electroweak sector of the minimal standard model with all 
known physical parameters, essentially G M , raw, m z and mtop, 
but parametrized with mjj- 



is of interest, fn cosmology, the neutrino degeneracy 
(the net neutrino number) of the universe is a poorly 
known number. Best limits are given by constraints from 
big bang nucleosynthesis and cosmic microwave back- 
ground radiation which limits the degeneracy parameter 
£„ = ji v /T v to £ MiT < 2.1 for the muon and tau neu- 
trinos and to £ e < 0.3 for the electron neutrino where 
/i„ are the neutrino chemical potentials and T v is the 
temperature of the neutrino background If such 

large chemical potentials were present in the very early 
universe then it raises a question about how they affect 
the electroweak thermodynamics and especially the elec- 
troweak phase transition. It has, for example, been pro- 
posed that the presence of a large lepton number asym- 
metry might explain the absence of topological defects 
h^i as well as the observed baryon number asymmetry 
[l3j. Finally, a comparison between QCD thermodynam- 
ics and electroweak thermodynamics is interesting. QCD 
thermodynamics has of course attracted a lot of interest 
during the last years due to experiments carried out at 
the moment at RHIC in Brookhaven and in the future 
at LHC in CERN. It should be interesting to see how 
the properties of the QCD phase transition as a function 
of baryonic chemical potential and number of light fla- 
vors (strange quark mass) compare to the properties of 
the electroweak phase transition as a function of leptonic 
chemical potentials and number of light bosonic degrees 
of freedom (Higgs mass) . 

The role of finite lepton number density in the thermo- 
dynamics of the electroweak theory has been discussed al- 
ready in the literature HUH E BUj. Those 
studies rely on perturbative one-loop calculations of the 
effective potential and the conclusion made is that the 
critical temperature increases with increasing chemical 
potentials. This can be understood in terms of Bose- 
Einstein condensation of the Higgs field due to finite 
chemical potentials related to gauge charges. The fate of 
the W ± boson condensate, predicted in [13, at high tem- 
peratures is also discussed on the same footin g Il8| . 
Vector boson condensation is also discussed in [21| 

Purely perturbative calculations are, however, doomed 
due to infrared divergences and a nonperturbative study 
is needed, in general. However, direct Monte Carlo stud- 
ies of the full electroweak theory at high temperatures 
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and finite chemical potentials are very difficult for nu- 
merous reasons. For example, the system is character- 
ized by a multitude of scales extending from irT (mass 
scale of nonzero Matsubara modes) to g 2 T (mass scale of 
the magnetic sector of the system). This leads to a need 
of large lattices in solving the properties of the system. 
Furthermore, chiral fermions are notoriously very hard 
to implement on a lattice. At finite density there is also 
the famous "sign" problem: the fermionic determinant is 
complex and thus the integration measure is not positive 
definite which spoils importance sampling. 

We approach this problem by generalizing the success- 
ful methods of 0, 0] to finite chemical potentials. That 
is, we calculate the dimensional reduction of the full 
4-dimensional electroweak theory at high temperatures 
and finite densities to a 3-dimensional effective field the- 
ory. The effective theory can then be solved using Monte 
Carlo methods in order to find out the phase diagram. 
This method is, however, by construction limited to small 
chemical potentials. 

The paper is organized as follows. In section^ we de- 
fine the theories and give the matching between them. In 
section ITTT1 we give the results for the phase diagram and 
in the final section, section llVl we discuss their meaning. 
Results for the required sum-integrals are given in the 
Appendix [X] 



II. DIMENSIONAL REDUCTION OF THE 
ELECTROWEAK THEORY AT FINITE 
CHEMICAL POTENTIALS 



is taken to have a nonzero Yukawa coupling. The conven- 
tion for the Euclidean gamma matrices is as given in [|| . 
The bosonic fields (<fi) are periodic in t while fermionic 
fields (t/0 are anti-periodic. Thus they can be expanded 
in Fourier series (Matsubara modes) 

oo 

4>(t,x) = M*y 2n7rTT , (2) 

n— — oo 
oo 

V>(t,x) = Mx> l(2n+1)7rTT - (3) 
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We will employ the power-counting rules g' 2 
A ~ <7y ~ g 2 . The Lagrangian is CP symmetric, 
culations are performed in Landau gauge. 

As well known, the thermodynamics of any system is 
described by the partition function defined as trace of the 
density matrix 



Z = Tr e -^ H -^ N » 



(4) 



Here Nk are all the conserved (global or local) charges 
of the system and Hk are the corresponding chemical po- 
tentials. In the electroweak theory, at classical level, the 
lepton number currents and the baryon number current 
are conserved independently. However, due to the trian- 
gle anomaly, these currents are not conserved in quantum 
theory (for a review see, e.g. j^) 

d^j^ oc g 2 e a p^G a al3 G a ^ for each current. (5) 



In this section we describe the construction of high 
temperature effective field theories at finite chemical po- 
tential. 



It is thus possible to form only n f conserved linear com- 
binations of these currents. These are usually defined to 
be 



X, 



A. The fundamental theory 

The electroweak theory at finite temperatures is de- 
fined by the Euclidean action 







with 



S = J dr Jd A xC 

C = (D^D^-v 2 &$ + \(&$) 2 + ~G a ^Gl v 

^F^F^+IlIPIl + e R pe R + q L flq L (1) 
+u R Ipu R + d R pd R + g Y (q L $t R + i R &q L 

Here D M = d^ + IigA^T 11 + Yig' where / and Y are the 
weak isospin and weak hypercharge of the corresponding 
doublet /singlet, II and q L denote the left handed lepton 
and quark doublets and e R , u R and d R denote the right 
handed leptons, up type quarks and down type quarks, 
respectively. Also, = d^A* - <9„A£ - gt ahc A\Al, 

Ffj, u = d^By - dyB^ and $ = ir 2 $*. Only the top quark 



-B — Li, i = 1 ... n/ 



(6) 



where n/ is the number of families, B is the baryon num- 
ber and Li are the lepton numbers for each family 



B 

U 
where 



gX) j d3x 5c,/7o9c,/, 

d 3 x (ei7 ei + Vi^a L Vi) 
a L = ^(1-75)- 



(7) 



Here / and c stand for flavor and color, respectively, and 
q c .f are the quark fields. The remaining current njB + 
J2iF% is not conserved. 

In addition to these globally conserved charges there 
are locally conserved charges related to the gauge sym- 
metries of the theory. Of the four gauge generators we 
can choose two mutually commuting ones for which it 
is possible to assign chemical potentials. One must be 
the hypercharge and as the other one it is convenient to 
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choose the third component of the isospin. The corre- 
sponding currents are 

■y lv-/l- 4_ 2 T , 

fam. 



fam 
2 



- ((£> M *)tr 3 * - $V 3 D M $) - e 3bc A^ b G c ^. (8) 

The sums run over the families. Chemical potentials re- 
lated to gauge charges cannot, however, be chosen freely. 
In thermal equilibrium the system must be neutral with 
respect to gauge charges. This requirement fixes the val- 
ues of these chemical potentials, which are then functions 
of temperature and chemical potentials related to global 
charges. This can be seen explicitly below. 

Taking all the conserved currents into account, the par- 
tition function is given by the path integral |18| 

Z = Trexp(-/3(i7-/i,X l -/iyQy- MT 3Q T 3))(9) 
-is- Tdr^^X^ (10) 



I Vip exp 
= / Dip exp 



,/3 «£ y 

J ° i=i /_ 

,/3 / »/ 

-S + / dr fi B B + y2 VLiLi 



where ip denotes the set of all the fields. Here we have 
defined (<p excludes B and Aq) 



/'B 



rif ^-^ 

1 i=i 



S = S 



<P,B 



WY A3 



(12) 



with S containing also the gauge fixing and ghost terms. 
We see that the constraint n//us + J2i = is satis- 
fied. We also note from Eqs. (jTU|) and (|T2*)l that after the 
integration over the conjugated momentum fields is done 
(going from Eq. © to Eq. I|10| l). the chemical potentials 
related to the gauge charges enter the path integral the 
same way as the static modes of the temporal components 
of the corresponding gauge fields. Thus, we can inter- 
pret them as acting as backgrounds for the gauge fields. 
Therefore, writing B -> B + (B ), A 3 -> A% + (A 3 ,) 
and requiring stationarity of free energy with respect to 
expectation values (Bq) and (Aq) (condition for thermal 
equilibrium) is equivalent to requiring neutrality with re- 
spect to corresponding gauge charges: 







<91nZ 

d(Bo) 



,d\nZ ig'Q\ 



T 



(13) 



and similarly for (Aq). It is then convenient to redefine 
B and A\ in Eqs. (JTDJ and JT^J) as 



This way the action in Eq. (|12(l becomes the standard 
electroweak action for these redefined fields and fiy and 
fix 3 are not explicit in the path integral anymore. The 
chemical potentials can be recovered as expectation val- 
ues for the redefined B and Aq at equilibrium 2 



/n\- lfiY /A-^ 3 
~ ~ : ~ ~T~' 



(15) 



This explicitly shows that these chemical potentials can- 
not be chosen freely, but are fixed by requiring that the 
system be in thermal equilibrium. 



B. Dimensional reduction at finite fi 

Due to infrared divergences which arise when integrat- 
ing over the bosonic zero modes 4>o(x) (static modes 
in r), the path integral (|10fl cannot be reliably evalu- 
ated within perturbation theory. The reason for this is 
that these modes are light when the temperature is much 
larger than any other mass scale in the theory and there- 
fore the high temperature expansion parameter g 2 T/E 
is large for them. All the other modes (nonstatic in r), 
4> n ^o(x) and ip n (x), are on the other hand always very 
massive, m ~ irT and can therefore be integrated out per- 
turbatively, as can also the static modes with \p\ > gT. 
We are then led to a very natural idea of formulating 
a three dimensional effective field theory for the static 
modes <f>o(x) 0113]. This effective theory is defined to 
be the most general theory for the static modes respect- 
ing the required symmetries. It reproduces the Green's 
functions for these static modes to a controllable accu- 



and (11) racy. 



Finite fermion number density affects dimensional re- 
duction in two ways. First, the renormalization of the 
fields and parameters as the heavy modes are integrated 
out changes when compared to the case [i = 0. Second, 
the symmetries of the fundamental four dimensional the- 
ory are reduced which gives rise to terms in the effective 
theories absent at /i = 0. More precisely, the introduc- 
tion of chemical potentials to the theory leads to terms 
in the path integral which break C but preserve P and 
T thus making the theory, in addition to being C and P 
breaking, also CP and CPT breaking. The effective the- 
ories may therefore contain terms which break CP and 
CPT and which do not appear at fj, = 0. Such terms 
must nevertheless still preserve 3-dimensional gauge- and 
rotational invariance as well as T invariance. 

The first effective theory is obtained after integrat- 
ing out the nonzero Matsubara modes. The resulting 
effective theory is a 3-dimensional SU(2) x U(l) gauge 
field theory with a fundamental scalar doublet (Higgs) 



g'B Q -> g'B - i/j, Y , gA% 



9Al 



(14) 



2 To be precise, this requires that the expectation values of the 
original Bo and Aq vanish in equilibrium. This certainly is the 
case when the system is neutral. 
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FIG. 1: Diagrams leading to new terms in the effective theo- 
ries. Solid lines correspond to fermions, wavy lines to gauge 
fields and dashed lines to scalars. 

and four adjoint scalars corresponding to the tempo- 
ral components of the gauge fields in the fundamen- 
tal 4-dimensional theory. The dimensionally lowest or- 
der CP and CPT violating terms arising from the finite 
chemical potentials in this theory are (in 3-dimcnsions 

[<P] = GeV 1 / 2 ) 

dim = GeV 5 : iBq, 

dim = GeV^ : iB%, i& A%T a $ , i&B <5>, iB Q A%A%, 

dim = GeV 2 : <,,,,!!, F jk , 

e ijk (AfG] k - '-ge^AfA^Alj . (16) 

which arise from diagrams given in Fig. ^ The factors 
of i are chosen in such a way that these operators are T 
invariant (with T transformations adopted from the 4- 
dimensional theory). The effective theory does not con- 
tain the terms TtAqT - and Tr(AgT a ) 3 since these vanish 
identically due to the properties of the SU(2) generators. 
The coefficient of the term _Bq is of the order g' 3 ~ g 9 ^ 2 
and when the matching of the Green's functions is done 
to order g it can be neglected. 

The possibility that there are terms linear in the ad- 



joint scalars in the action of the effective theory is note- 
worthy. Such terms induce condensates of the corre- 
sponding fields in equilibrium. As already discussed, 
such condensates for the adjoint scalars are equivalent to 
nonzero chemical potentials for the gauge charges. Thus 
the emergence of linear terms in the effective theory takes 
care of neutrality of the system with respect to gauge 
charges. 

The last two terms of Eq. I|16f) . the so called Chern- 
Simons terms, are interesting. There of course is a vast 
literature on the physics induced by them (for a review, 
see [24]]). That such terms appear in chiral gauge field 
theories when fermions are integrated out was first ob- 
served in |25|. However, in the present study we ob- 
serve that the coefficients of these terms in the effec- 
tive theories vanish due to the nonconservation of B + L 
(n//z_B + Yli Mi< = 0)- Thus they do not play any role in 
them. The role of those terms at smaller temperatures, 
where they may be important, has been discussed in |26| . 

The second effective theory is obtained after further in- 
tegrating out the adjoint scalars, the zero components of 
the gauge fields. The resulting theory is a 3-dimensional 
SU(2) x U(l)+Higgs gauge field theory. The form of the 
theory is fully determined by the gauge invariance. There 
cannot, therefore, be any terms in the theory that would 
be absent at p = as far as symmetries are considered 
(apart from the Chern-Simons terms, which nevertheless, 
as stated above, turn out to be absent). Finite chemical 
potentials show up only in the mapping of the parameters 
of this theory to the physical variables. 



C. Integration over the superheavy modes 

The first effective theory in its most general form is 
defined by the Lagrangian 



£1 = -G%G% + -FijFij + (A* ) f A* + m 2 $t$ + A 3 (&$) 2 + - (A^S) 2 + 2 m D A o A o 

+ ±\ A A a A a A b A b + \ (d.Bof + \m%Bl + h 3 &$A a A% + h^$B% + ~g 3 g 3 B &A%T a $ (17) 
+ae ijk yAfG" k - l -g a e ahc A1A)Al^ + a'e m BiF jk + Kl B + p&A a r a <Z> + p'&$B + PG B A a A a Q , 

I 



where A^o = d l A? ) -g 3 e abc A l >A'jj. The parameters of this 
theory are to be matched to those of the 4-dimensional 
theory up to order g A . The factors of i are included in 
the coefficients n\, p, p' and pq. 

In general, theories are related by matching corre- 
sponding Green's functions calculated in each theory. If 



the fields are renormalized as (cp denotes a generic field) 

fin = ifZ^l-D, (18) 
then the TV-point Green's functions are related by 
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given in Eqs. (141), (142) and (143) of we get 



FIG. 2: The required diagrams for Debye masses and the 
couplings. Solid lines correspond to fermions, dashed lines to 
scalars and wavy lines to gauge bosons. 




V 3 / 

^(f) + E-^) 



i=l 



(21) 



O O 



FIG. 3: The required scalar 2 point functions for the mass 
parameter. Solid lines correspond to fermions, dashed lines 
to scalars, wavy lines to electroweak gauge bosons and curly 
lines to gluons. Crosses denote counterterms. 



Here we have defined the function A{\i) (see Appendix 



<m = v> ( 5 



2^f 



2^T 



2j E + 2 In 4. 

It is now straightforward to find the modifications of the 
renormalization of the parameters. The Debye masses 
get a correction coming from the diagrams in Fig. J5| 
These are needed only to order g 2 (as discussed in (6j) 
and thus the modification of the field renormalizations 
do not affect the Debye masses. Denoting by subscript 
fx = the results from Eqs. (160) and (161) of we get 



where i labels the different fields and is the number 
of times the field i occurs in the Green's function. This 
matching has been performed for the minimal standard 
model at zero chemical potentials in [|| . At finite fermion 
number density those results are modified. 

Let us denote by Ar^ and AZ V the change in the 4- 
dimensional Green's functions and field renormalizations 
due to finite chemical potentials, and by Ar 3 rj the change 
in the 3-dimensional Green's functions due to changes in 
the parameters of the effective theory. We then get from 
Eq. JED that 



Ar ({« 

^ L 3D 



}) _ rpN/2 



n 



o.-nj/2 \ Ar ({n 
- S 'i,Ai=0 I ^ L 4D 



E 



az, r 



4D 



(20) 



which holds when Z v are calculated up to one-loop or- 
der. This formula allows us to calculate the changes in 
the mapping of the parameters of the effective theory to 
physical variables. 



1. Changes in the coefficients of the terms present already 
at \i = 

The field renormalizations can be calculated from the 
momentum dependent part of the two-point functions. 
Denoting by Z° the results at zero chemical potential 



' D 



m' 2 




(23) 



The coupling constants, on the other hand, are modified 
due to both changes in field renormalizations and changes 
in loop integrals. The required diagrams are shown in 
Fig. [21 The results are 



A, 



A 



Sw-^M(f) 



9t 



#3,^=0 + 




"(f) 




where A 3 ^ =0 , h 3..p,=n_ and g 2 yfJi=0 are given in Eqs. (150), 
(147) and (146) of ||, respectively. One may note that 
Xa, the self coupling of the adjoint scalars, does not get 
any corrections from the chemical potentials and thus Eq. 
(162) of holds. Also, h' 3 = g'^/A. 

Last, the scalar mass parameter is calculated to two 
loops. One must carefully obtain all the contributions 
from fields renormalizations and loop integrals. The di- 
agrams needing recalculation are shown in Fig. The 



result is m|(/x) = 
given by Eq. (156) of 



o(M) ' 
and 



Am 2 where m 3 „. 



:o(m) is 
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Am 



12tt 2 



Ms 



3.9 



16tt 2 
9 



AT 2 _ S^T 2 _ 5yT^\ fff-Ml r 

2 16 4 J V 3 / 64tt 4 



<7 2 L 6 (M)-<£ (L f {fi)-A(ff)) 



9<4 + + 16^ 5s 2 ) -4 (^) - + \ 9 W ~ ~ Wydl ~ IBS 



? 4 1 ^ 2 2 io\2 1 c 2 2 27 4 
y + j9y9 - 18Ag y - lGg Y g s ~ ~^9 



9g Y g 2 - 6 (3<& - 8^ 5s 2 ) L b (/2) + 9<£Z/(M) 



48.3 -"Vr(| + il 

9 



2 1 ' IGgWs) (41n2-l) 



~ 2 9y9 



27 



E 



8tt 2 



#(mlJ + -T7T-3" ln 



/ 1 1 i^Lj 
\2 ' 2ttT 




T 2 



1287T 2 



Mb 



11527T 4 



(25) 



r 



Here the function 23(/i,) (see Appendix lA"|l is defined by 



+ C '(_l i_ J^i. 
+ ^ 1 '2 2vrT 



(26) 



and the functions and Lf(p,) are as defined in 



Lft(jS) = In ^ - 21n47r + 2 7B 
u 2 

L f (p) = ln^-21n7r + 2 7£ . 



(27) 



Here p, is the renormalization scale in the MS scheme. 
Note that Am 2 is independent of M when the running of 
g\ is taken into account. 



The corresponding diagrams are given in Fig. The 
coefficients of the other new terms are needed only to 
one loop order. The result for them is 



9 
1 

9 

9G 

a 



_J_ 5 y T 3/2 ( MB 



9 



8tt* 
2 



ttT 

n f 



32tt 2 

J2 



9 ' 



32tt 



n /MB +J2^ L 'j = °' 

V i=l / 

■ n f9B + E^ L » = °" 



(29) 



2. ATew terms 

As already pointed out, chemical potentials also induce 
new terms to the effective theories. These arise from the 
diagrams in Fig^ The most important one of these is 
the term linear in Bq which is related to the neutrality 
of the system. Calculating to two loops (not including 
contributions of the order g' 3 ~ g 9 ^ 2 ) we get 



[ 9 / T 5/2 



9ff 2 \ y> Ml, 
64tt 2 I ^ ttT 

' 1=1 



- 1- 



5g| 
32tt 2 



9g 2 
64tt 2 



9b 
2tt 2 / ttT 



1 / Mb\ 2 
9 VttT/ 



These come from the diagrams in Fig. As already 
noted, the Chern-Simons terms vanish due to the non- 
conservation of B + L which sets nffis + 53« ML; = 0- 



(28) 



FIG. 4: Diagrams contributing to the linear term. Solid 
lines are fermions, dashed lines scalars, wavy lines electroweak 
gauge bosons, curly lines are gluons and external legs are -Bo 
legs. 
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D. One-loop effective potential 

Before integrating over the adjoint scalars A$ and 
Bq it is instructive to consider the effective poten- 
tial V cS (ip, (Af t ), (B )) for the condensates ((<&) = 
l/y/2 (0, y) T ). Although not completely reliable, V e e is 
known to give a rather good description of the phase tran- 
sition at small ran where the scalar self coupling is small. 
Studies of the electroweak phase transition at finite chem- 
ical potentials using perturbatively derived effective po- 
tentials have been carried also previously [Til IT7L ITsj . 

At this point we are only interested in the 
qualitative effects finite chemical potentials can 
have. Therefore, in order to simplify the proce- 
dure, we neglect the contributions from the terms 
X A A^A b A b , p$Ug r Q $, p'B $t$ and PG B A a A a . 



This is motivated since ~ g 4 and pq ~ g 7 / 2 are 
small and since the three point vertices Bo&& ~ p' 
and ~ p are negligible when compared to similar 

vertices obtained from the four point vertices Bq&Q 
and Bq& AqT 11 ^ after annihilating a Bo leg by the K\ 
vertex. Likewise, we will use /13 = <?|/4, h 3 = g 3 2 /4, 
neglecting higher order corrections, which simplifies 
some expressions below. This is also consistent with all 
the approximations above. Furthermore, we will only 
consider quantum fluctuations of the magnetic sector 
Af, Bi and of the scalars Ag, Ag, and treat all the 
condensing scalars Bq and Ag only at the tree level. 
This approximation is adequate to show the effects 
of chemical potentials. A standard calculation of the 
one-loop effective potential gives 



V cff (^ (Bo), (AD) = -m 2 ^ 2 + -V + ^l(Al) 2 + -mg(So) 2 + 5 (h 3 (A%) 2 + h' 3 (B ) 2 )^ 2 - -g 3 g' 3 ^(B )(Al) 

+*i(Bo) (jffaV + 4<^) 2 ) 3 / 2 + 3 -(g 2 + g' 2 f' 2 ^ + 2(ml + h^f'^ . (30) 

I 



Requiring neutrality with respect to gauge charges en- 
forces the conditions 



dV c g 
d(B ) 







dV cS 
d(Al) 



0. 



(31) 



which gives us 

(Bo) 
(B ) 



Bo)o + {B ) u (Al) = (A 3 ) + {A 3 )i with 

, . h?. \ 1 
,- I 1 + - 

TO, 



n 2 

-TV 

L D 



(Al)o 
(Boh 



(Al) 



Kigzg'zv 2 



(hi , h' 3 \ 



3 gig' 3 

16-7T m 2 D m'p 



(AlW 



v^njWo 



(32) 



"D 



Here (B ) and (Aq)o are tree level contributions to the 
Bo and Ag condensates, and (-Bo)i and (Aq)i are cor- 



rections to those from the one-loop term of the effective 
potential. Inserting these back to the effective potential 
in Eq. H3U|) gives us then the effective potential for the 
Higgs expectation value 



2 "" 
1 

+ 4 



3 - 7T h $ mD + 

Z7T 



h! 3 n\ 



1 



hi 



2 hr^ t\i 



47T TO£) 



h 
m 



r 

m, 



32tt 



2gl + (g 2 +g' 2 f' 2 )^ + 0{^) 



(33) 



The error made by neglecting higher order terms in ip is 
small near the phase transition as can be seen from the 
example in Fig. |3J In terms of the physical parameters we 
get (« corresponds to only tree level matching between 
couplings of the 4-dimensional theory and physical vari- 
ables, given for example in Eqs. (184)-(185) of 6]. The 
exact one-loop matching, given in Eqs. (183) and (194) 
of 6], would not change the qualitative considerations.) 
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VJ<p), H = 30 GeV, T = 99.0 GeV 




(p/GeV 

FIG. 6: The evolution of the effective potential at fixed tem- 
FIG. 5: The one-loop effective potentials at the phase transi- perature T = 98 GeV and Higgs mass m H = 60 GeV as the 
tion point for m H = 60 GeV. As can be seen, omission of the chemical potential is increased, 
higher order terms (see Eq. (JSHJ) has only a small effect. 



where we have set all the leptonic chemical potentials 
equal and fi = (il, = —^b and only leading order terms 
in fx are kept. 

This effective potential gives a qualitative picture of 
the effect of the chemical potentials. First, the critical 
temperature increases due to decrease of the Higgs mass 
term. Second, the scalar self-coupling increases leading 
to a smaller Atp at the transition and thus to a weaker 
transition. Also the "barrier" responsible for the first 
order phase transition is lower. All these effects can be 
seen in Fig. [5] In Fig. [S]it can be explicitly seen that 
finite chemical potentials tend to break the symmetry of 
the theory. 

The leading /i-induced corrections to both the scalar 
mass parameter and self-coupling come from the Bq con- 
densate. Thus, at leading order the leptonic chemical 
potentials change the properties of the phase transition 
through generating nonzero chemical potentials for the 
gauge charges. These couple to the Higgs field and thus 
change the dynamics of the Higgs field. 

We can also note that the W ± boson mass is reduced 
in the broken phase due to the A$ condensate given in 
Eq. (Eg, 3 

2 1 / 2 576 p 2 2 \ 2 
mw * 4 [ 9 * " 146412^ (35) 

If fi/T is large enough, the bosons may become un- 
stable leading to a M /± -condensate. We, however, restrict 



3 This means that the physics behind the condensation is 

related to the nonzero chemical potentials for the gauge charges. 
These chemical potentials couple to the W^ 1 bosons and this 
condensation is nothing but Bose-Einstein condensation due to 
these chemical potentials. 



ourselves to small chemical potentials and study phenom- 
ena near the phase transition where ip 2 <C T and then the 
condensation is not relevant. 

Although the perturbatively derived effective poten- 
tials serve well in giving a qualitative, and even a quan- 
titative, picture of the phase transition, they cannot be 
trusted as the Higgs mass increases. This is clearly seen, 
for example, in predictions of the nature of the phase 
transition. Perturbative calculations predict a first order 
phase transition for all Higgs masses while nonperturba- 
tive studies show that there, in fact, is only a crossover 
for run > 72 GeV at /i = 0. Since the physical Higgs 
mass is large, it is important not to rely on perturbative 
calculations. 



E. Integration over the adjoint scalars 

The phase transition occurs when the scalar mass pa- 
rameter m 2 becomes small, m 2 ~ g^T 2 . Then the adjoint 
scalars may be considered heavy, m 2 D ~ g 2 T 2 , and they 
can therefore be integrated out. The resulting effective 
theory is defined by the Lagrangian 

£2 = ^Gtj+^FijFij+iDi^ A$+m 2 $t$+A 3 ($t $ ) 2 , 

(36) 

This reduction step C\ — > £2 is performed in jfj in the 
case Ki = p = p' = pa = 0. We now generalize this 
to finite K±, p, p' and pq. We use as the expansion 
parameter /13 jm d and the goal is to calculate the correc- 
tions up to order (h^/mD) 2 "m 2 D for the scalar mass and 
to order (/i3/mi)) 2 mz) for the couplings in the reduction 
Ci — > £2- To keep track of the contribution of different 
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(a) 



(I.) 



FIG. 7: (a) A tree level diagram in the theory given by L\ 
leading to higher order term (&&) 4 in the last effective the- 
ory, (b) Leading order contribution to the Green's function 
( (<3>^ *E > ) 4 ) in the theory given by €2- (c) A tree level dia- 
gram in the theory given by £1 leading to higher order term 
<E>) 2 A% Af in the last effective theory responsible for the 
W condensation . Dashed lines correspond to the fundamen- 
tal scalar, solid lines to adjoint scalars, wavy lines to gauge 
fields and dots to the Ki vertex. 



terms we use the following power counting rules 

2 7 ^3 12 U ^ 3 u ( ^ 3 , 

g ~ h 3 ~ m D , g 3 ~ h 3 ~ h 3 ~ m D , 

12 ^3 2 ( ^3 \ a 5/2 

m D m D , Ki ~ m£ , 

m D \ m D J 

( h 3 \ Q_1 / 2 3 , 2 / ft 3 \" 3/2 

/ /i 3 \ 2+ " 3/2 
\m D J 



(37) 



These arise from the power counting rules of the original 
theory supplemented by setting /z/(ttT) ~ (h 3 /rriD) a for 
some a where p can be any of the chemical potentials. 
Due to large mass of the top quark we have also relaxed 
the power counting of p and p' by treating g\ ~ 1. This 
protects us from neglecting p and p' in situations where 
they would be important. 

It is essential to keep the chemical potentials suffi- 
ciently small in order to keep this last reduction step 
meaningful. Consider, for example, the set of diagrams 
in Fig. [7^. Such diagrams would lead to terms of the 
form (<J>t$)", n > 3 in the final effective theory. If such 
higher order terms are not to be included in the effec- 



tive theory, then the contribution of those terms to the 
corresponding Green's functions must be sufficiently sup- 
pressed so that they can be neglected. More specifically, 
in the theory described by £2, the leading contribution 
to the Green's function ((<I>^$)") comes from the scalar 
loop in Fig. [7jp and is by naive dimensional calculation of 
the order \ 3 m\~ 2n ~ (h 3 / 'm D ) 3 ~ n m 3 ] ^ n near the phase 
transition where rh 3 ~ {h 3 /m£ ) ) 2 'm 2 D . The contribution 
from the graph in Fig. [TJl on the other hand is of the 
order of (h 3 /m D ) 2a+n ~ 3 m 3 p n . Thus, if we require that 
the effective theory can reproduce the Green's functions 
up to order g 3 ~ (h 3 /mn) 3 and we neglect the higher 
order operators produced by the graphs in Fig. [TJi, then 
we must require that 



. 3 -in M / 

a > — or equivalently — — 
2 7rT 



ff. (38) 

This determines the powercounting rules for m, p, p' 
and pa which are used when integrating out the adjoint 
scalars. 

Another set of interesting tree level diagrams are those 
in Fig. 0:. They would lead to terms of the form 
A? in the effective theory £ 2 - These are in- 
teresting since it is these terms that are responsible for 
the condensation, as can, for example, be seen from 
Eq. H35fl . There the leading correction to the mass 
is ~ (p and thus the term in the effective theory £ 2 
that would be responsible for this correction would be 
~ (<E>i $) 2 A c [ Af . The above determined power counting 
rules allow us to neglect these terms but this means that 
the effective theory £ 2 cannot predict W condensation. 
As already discussed, this is not a problem when we study 
phenomena near the phase transition and at small chem- 
ical potentials. 

Since the reduction at k\ = p = p' = pc = is given 
in |(J we now only need to take into account the contri- 
bution from the new terms. Calculating to the accuracy 
mentioned before, they only contribute to the scalar mass 
parameter (up to one-loop level) and scalar self coupling 
(up to tree level) . Other contributions are of higher or- 
der. The required diagrams are shown in Fig. |SJ The 
results for the parameters are 



™2 



fh 2 



3.0 



+ 



h' 3 K 



D 



A 3,0 — 77 



2 p'm 


1 


m% 


All 


ah? 4 


Ah' 2 \ 


Am% 





Ami 



Ah' 2 



D 



D 



'1 



Ah' 3 p> 



rn 



D 



3g3g' 3 p . 4ti 3 p' 



m D 



111: 



12 



D 



«1 
m D 



m D 



(39) 



Here m§ and A^o are as given in Eqs. (174) (first equal- 
ity) and (169) of 6} (where they are denoted by mfand 
A3). All the other couplings of £2 are as given in [6j as 



functions of the parameters of C\. 

The construction of the theories C\ and C 2 differs in 
one qualitative aspect. Although neither theory is appli- 
cable at high chemical potentials, p > ttT, there never- 
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FIG. 8: Additional diagrams needed in integration over the 
adjoint scalars. The dashed lines correspond to the funda- 
mental scalar, solid lines to adjoint scalars and the dot is the 
Ki vertex. 



theless was no expansion in li/(ttT) in the construction 
of C\. That is, the matching of the parameters of L\ 
to those of the 4-dimensional theory is given to a cer- 
tain accuracy in g 2 for arbitrary \i. This is not true for 
the matching between C\ and £2- There it was essen- 
tial to assume that the chemical potential is small. This 
is easy to understand since some of the couplings of C\ 
arc directly proportional to n/(jrT) and thus /i must be 
small when C\ is studied perturbatively. We may there- 
fore assume that C\ is applicable to somewhat higher 
chemical potentials than £2, up to [i < ttT. In the range 
g'irT < ix < irT the dynamics of L\ may be dominated 
by nonperturbative effects. 



III. THE PHASE DIAGRAM 

The effective theory for the light modes is infrared di- 
vergent in perturbation theory and thus it cannot be re- 
liably studied perturbatively. It has, however, been stud- 
ied nonperturbatively by Monte Carlo studies in 0. 

The theory is parametrized by four parameters. It is 
convenient to express three of them in a dimensionless 
form while leaving one of them to give the energy scale. 
We define 



A3 

w v 



~4 



_~9? 



% (40) 
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and leave g 2 to give the dimensions. The value of z 
is essentially fixed by the Weinberg mixing angle, z ss 
tan 2 W ~ 0.3. The value of y is tuned to find the phase 
transition at a fixed value of x which determines the na- 
ture of the phase transition. 

The phase diagram of the effective theory is given in 
Fig. [5] The continuous line is a curve fitted to the lat- 
tice results which are given in [7J and 9] for SU(2)+Higgs 
gauge theory. The effect of the U(l) subgroup is to in- 
crease the critical y slightly Q • The critical line given by 
a perturbative calculation is given by the dashed line. As 
can be seen, the perturbative result gives quite a good 
estimate for the value of y at the transition, y = y c {x), 
for small x. It, however, fails completely in describing 
the nature of the transition at high x. Perturbation the- 
ory predicts a first order phase transition for all x, while 
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FIG. 9: The phase diagram of the effective theory. The results 
are from and 0. 
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FIG. 10: Behavior of x c and y c as a function of /j, for differ- 
ent Higgs masses. The regions above the line x c = x 
and below the line y c 
region of the theory. 



V, 



endpoint 



endpoint 



correspond to the crossover 



Monte Carlo studies have shown that there is a first order 
phase transition at small x but that the first order phase 
transition line has a second order endpoint at x w 0.0983, 
y ss —0.0173 and for larger x no phase transition is ob- 
served 0)0- Thus, there is no phase transition for suffi- 
ciently large Higgs masses. 

The phase diagram in Fig. EJcan be expressed in terms 
of the physical parameters by the mapping described in 
section [H] For simplicity we set all the leptonic chemical 
potentials to be equal to each other, fi^ i = ji = —fiB- 
The theory is specified by giving the physical param- 

10- 5 GeV" 2 , m w = 
: 91.19 GeV, m t = 174.3 GeV and 
The Higgs mass is left as a free pa- 



eters the values G M = 1.664 
80.42 GeV, m z 
a s {m z ) = 0.118. 
ramctcr. 

The exact relation between (mji, T, n) and (x, 



y) 
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is a complicated function of //. However, the essential 
features of the effect of the chemical potentials can be 
seen quite easily by just taking the leading corrections 
due to finite /i into account. We get 



0.18 



yip) 



x(0) 

y(o)- 



1 96 



1 96 



g 2 1331 T 2 
1 16 ^ 2 



ff 2 1331 T 2 : 
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121 T 2 



1 



2 5 4 T 2 

1 16 /x 2 
7l21 T 2 ' 




(41) 



Thus the effect of finite /i is to increase x and decrease y. 
Perturbatively, at tree level the phase transition occurs 
at y — y c = 0. The critical temperature can then be 
solved to be 



Tn = 



1 



8A + 3g 2 + g' 2 + 4e£ 



8m 



H 



256 ,> . 
121 ) (42) 



which explicitly shows that a finite chemical potential 
increases the critical temperature. At one-loop order the 
critical line is given by y c x = l/(1287r ) and expanding 
this around the tree level solution T c = Tq we get for the 
critical temperature (to first order in (T c — T )/T ) 



T r = T n 1 



1 



16tt 2 8A + 3g 2 



4. An^ , 96 uf 



(43) 

The factor multiplying the tree level critical temperature 
decreases as the chemical potentials are increased but 
that decrease is negligible in comparison to the simul- 
taneous increase in To. Thus the critical temperature 
increases also at one-loop order. 

The above reasoning gives a valid qualitative under- 
standing of the behaviour of the system but to obtain 
quantitatively more reliable results we must map the 
phase diagram in Fig. [!5]to (T, /i, run) using the com- 
plete results from section [H] Then the values of x and 
y along the phase transition line, (T c ,fi,m H ) and 

y c = y(T c , n, mn), are given in Fig. E3as functions of 
/j. It can be observed that as the chemical potentials 
are increased, the subsequent increase and decrease of x c 
and y c , respectively, are fastest at small Higgs masses. 
This is easy to understand since x c and y c are essentially 
functions of n/T. Thus, at small Higgs masses where 
the critical temperature is lower, increasing [i leads to a 
larger increase in /i/T c than at large Higgs masses where 
the critical temperature is higher. Therefore the changes 
in x c and y c are also larger at smaller Higgs masses. 

This has an interesting consequence. The x c (/i) curves 
for different uih intersect and for sufficiently large chemi- 
cal potentials the value of x c is, in fact, a decreasing func- 
tion of the Higgs mass (at least for sufficiently small Higgs 
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(1 = 45 GeV 
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FIG. 11: The behavior of function of ran ■ 

masses). This can be seen explicitly in Fig. 1111 There- 
fore, under these specific conditions, the phase transition 
appears to become stronger as the Higgs mass is increased 
(again, at least as long as the Higgs mass remains suf- 
ficiently small). This, however, does not seem to have 
any physical relevance. Dimensional reduction is not re- 
liable at small Higgs masses where this effect is strongest. 
As the Higgs mass is increased also the chemical poten- 
tials must be increased in order to recover this anoma- 
lous behaviour. However, at these larger Higgs masses 
and chemical potentials, the value of x c is above the end- 
point value ^endpomt _ q.0983 and there is therefore no 
phase transition. Especially, at physical Higgs masses, 
m H > 115 GeV, there is no phase transition. 

The phase diagram in terms of the physical parameters 
is given in Figs. H2landlT3l The qualitative picture based 
on perturbation theory can be seen to be correct. The 
critical temperature increases with fj,. Some interesting 
thermodynamics can be deduced from this. According to 
the Clausius-Clapeyron relations 



dT 
dfx 



n b 



(44) 



where dT/dfi is measured along the phase transition line, 
and n Si b and s Sj b are the lepton number and entropy den- 
sities in each phase, respectively. Since dT/dfi is positive, 
the lepton number difference between the phases is of the 
opposite sign than the entropy difference. Furthermore, 
since entropy is higher in the symmetric phase (the high 
temperature phase) it means that the lepton number den- 
sity is higher in the broken symmetry phase. 

The behaviour of the endpoint of the first order phase 
transition line is also of interest. For a fixed Higgs mass, 
a finite chemical potential strengthens the scalar self- 
coupling at the transition point as seen from Fig. 1101 
Thus, increasing fi weakens the transition and the value 
of the critical Higgs mass where the first order phase tran- 
sition line ends in a second order endpoint is a decreasing 
function of /i. Especially, there cannot be a first order 
electroweak phase transition at physical Higgs masses in 
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FIG. 13: The electroweak phase diagram on the fi — T plane. 



the minimal standard model even when /i 0. Moreover, 
for sufficiently high chemical potential there appears to 
be no phase transition for any value of mn ■ The location 
of the second order endpoint in terms of the Higgs mass 
and critical temperature are given in table ^ for some 
values of the chemical potential. 



M 


end point 


rj-icnd point 
-* c 


GeV 


72 GeV 


109 GeV 


15 GeV 


71 GeV 


108 GeV 


30 GeV 


66 GeV 


104 GeV 


45 GeV 


52 GeV 


94 GeV 



TABLE I: The location of the endpoint of the first order phase 
transition line. 



FIG. 14: The phase diagram. Note that the temperature 
grows downwards 



IV. DISCUSSION 

In this paper we have determined quantitatively how 
the equilibrium phase diagram of the electroweak the- 
ory, parametrized by the Higgs mass, depends on the 
temperature and small leptonic chemical potentials. In 
particular, we have studied the change of the second or- 
der endpoint of the first order phase transition line. It 
is seen that the critical temperature increases and the 
critical value of Higgs mass where the first order phase 
transition line ends decreases as the chemical potentials 
are increased. These results are summarized in Fig. [Til 

It is interesting to qualitatively compare these re- 
sults for the critical temperature T — T c (r7ijj,/i) of the 
electroweak phase transition to the critical temperature 
T = T c (m s , A*b), m-s — strange quark mass, /ib = bary- 
onic chemical potential, of the QCD phase transition. In 
QCD there is also a first order critical line ending in a 
2nd order point: if \ib — then for m s — (Nf = 3) 
there is a first order phase transition while for m s = oo 
(Nf = 2) there is only a crossover. Thus, the situation 
on the m s — T plane is like in Fig. 1121 However, in con- 
trast to Fig. ED on the ^b~T plane the critical lines for 
different m s bend downwards starting from (is =0 (in 
Fig. I12l the /is > curves would reside below the /is = 
curve). As discussed after Eq. I|44|) . this implies that in 
the QCD phase transition the entropy density and the 
baryon number density change in a similar way. Also, 
in QCD, when m s is large enough, there is a first order 
phase transition only if fj,s is greater than some critical 
value /iB jC (m s ), again in contrast to Fig. 1131 which says 
that there is a first order electroweak phase transition 
only if n is less then some critical value /x c (mjy). 
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The effective theory used in this paper to study the 
electroweak phase transition cannot be used at large 
chemical potentials. Therefore, it cannot, for example, 
predict the emergence (or perhaps lack of it) of the 
condensate. The reason for this is that the driving force 
behind the condensation, Bose-Einstein condensa- 
tion, would appear in this theory through higher order 
terms of the form (&&) n AfAf, n > 2. Those terms 
would modify the mass in such a way that the con- 
densation might occur. We have, however, neglected such 
higher order terms from the theory which is possible if 
the chemical potentials are small. The condensa- 
tion is in any case a high \x phenomenon and thus our 
approximation is self consistent. 
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APPENDIX A: SUM-INTEGRALS AT FINITE fx 



sional momentum integration and then performing the 
sum using Riemann zeta functions. Denoting the sum- 
integrals by 



e^ £ V 



/ n— — oo 

where ft is the scale in the MS scheme, we get 



Po = (2n+l)irT—ifi, 
(Al) 



pi 



(P 2 ) a 



x2Re C(2a-2n-3 + 2e,-- ''" 



2 2ttT 



I; 



2n+l _ 



„2n+l 

Po 

(P 2 ) a 

r(a-3/2 + e) 
87r 3 / 2 r(a) 



x2Tm 



(e 7 V) e T(2^T) 



2n+4-2e-2ct 



C [2a - 2n - 4 + 2e, 



1 i/i 



2 2vrT 



(A2) 



In this appendix we give the results for the required 
fermionic sum-integrals for arbitrary \x. Similar integrals 
also appear in the context of computing QCD quark num- 
ber susceptibilities [27j |. 

1. One-loop integrals 

At one-loop we have two types of integrals. Both of 
these are easily evaluated by first doing the 3 — 2e dimen- 



where 

is the generalized zeta function. Its properties are dis- 
cussed below. Especially the following special cases, ex- 
panded around e = 0, are needed: 



7° - -^/l + A^ + e 
h ~ 24 1 , 2 T 2 



I l 2 



_3_/? 

^2^2 



r - 

in 



-24C' 



1 



1 

L <2 



-^-) + 24C 



-1 



(L / (^)-2 7B + 2-41n 2) 
1 i[i \ 



' 2 2ttT J 



16tt 2 
+0(e), 
iirT 3 fi 
12 7rT 



i + L /W -2 7 ,-41n2-^(i-^ 



1 



TT 2 T 2 



0(e), 



8ttitT V h 



0(e% 



(5- 


- i(x )} 




2irT J 



(A4) 

(A5) 
(A6) 

(A7) 



where the primes in zeta functions denote derivatives The results are exact in \x. Terms of the order e are 
with respect to the first argument and tjj(z) = d z hxT(z). 
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needed only for the first sum-integral. 



Some properties of the required special 
functions 



2. Two- loop integrals 

Most of the needed two-loop sum-integrals can be re- 
duced to products of one-loop sum-integrals. Only the 
sum-integral corresponding to the "setting sun" diagram 
must be evaluated explicitly. In calculating it we follow 
the procedure outlined in |28j . 

The integral we are interested in is 



1 



1 



Il f (k) (A8) 



,p f ,k„p 2k2 (p + k ) 2 Jk b 

Po = (2n + 1)ttT — ifi, kg = 2rmrT, n, m € Z, 

where we have defined the integral n^(fc) as in pgj . 
By going to configuration space where the propagator 
is given by (at e = 0) 



A(po,r) 



-\Po\r 



(A9) 



this integral becomes 

Il'(fc) = TJ2 Jd 3 re* k r A( Po ,r)A(po + ko,r) 



{po} 
T 



16n- 



{Po} 



ik r 



i\ r Q— \Po+ko\r 



(A10) 



Here we have used the notation of |28| where X){ Po } 
means summation over fermionic Matsubara modes. 

The influence of the chemical potentials resides now in 
the sum 



{Po} 



, borfcok = o -|fc |r ( C0s(2 M r) M \ 

v sinh(27rTr) 2irT J 



ik r 



d d r 



cos(2^r) |fcp| \ 
sinh(27rTr) 2ttT ) 



(All) 



|feok 



and thus we get 

(A12) 

From here we can proceed as in |2y by first subtracting 
the T = 0, = part (which contains the ultraviolet 
divergences) and then evaluating the remaining integral. 
The final result is 



64tt 4 



i a 

- + 21n^2 -41n47T + 2 



16.3 r(i + ^) 



+ 0(e). (A13) 



We were unable to find whether this integral had been 
calculated at finite fi previously. 



The central function is the generalized zeta function 
((z, q), defined in Eq. i|A3|l for Re(z) > 1, and its deriva- 
tive with respect to z, Q'(z,q). From the integral repre- 
sentation 



1 



T(z) 



dtt 2 



e* - 1 ' 



(A14) 



valid at Re(z) > 1, Re(g) > 0, it is possible to analyti- 
cally continue this function to the double complex plane 
(z,q), z ^ 1. For us it is sufficient to consider only the 
half plane Re(g) > where we can write 



1 



E 



Bn(q) (-1)" 

M-q)t 
dtt*- 1 ^— . 



T(z) 



(A15) 



Here B n (q) are Bernoulli polynomials. This expression is 
well defined as long as z ^ 1 and Re(g) > (the limit of 
l/(T(z)(z + n- 1)) is well defined at z = 0, -1, -2,...) 
and thus it can be used to define the analytic continuation 
of ((z, q) to that region. For the derivative we get 



C(z,q) = 



1 



Bn(q) (-1)' 
r(z)^r(n+l) (z + n- 



E 



r(z) 



dtt^lnt- 



3 (l-9)t 



1 



(A16) 



i>(z)((z,q)- 



All the required special functions appearing in the inte- 
grals can be related to these functions. This is trivially 
true for the one- loop integrals by Eqs. (|A2|) . The loga- 
rithm of the gamma functions present in the sunset in- 
tegral Eq. i|A13(l can on the other hand be written as 
(x = /V(ttT)) 
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( A1? ) 

Thus the mathematical properties of the non-polynomial 
part in /i/T of all the required integrals are given by the 
generalized zeta function. 

The asymptotic behaviour of the integrals at the limits 
fi/T — > and T/fi — > is of interest. The small fi limit is 
a straightforward Taylor expansion and is not given here 
explicitly. The large [i limit for the required functions 
is most easily obtained from Stirling's approximation for 
the gamma function: 
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(A18) 
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Thus 
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lnx 2 + O(x ). (A19) 



The last expansion can be derived with the help of Eq. 
(|A17|I after noting that the zeta function satisfies the 
relation 



^C'(z, <?) = -C(* + 1,9)- *C'(« + 1, g) (A20) 



and that £(0, q) — 1/2 — q. Low temperature expansion 
of the integrals (|A4|) , (|A5I) and (|A13|) is now straightfor- 
ward. We get 
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+ 0(e 2 ) 



(A21) 



Thus, although separate terms of the integrals diverge 
logarithmically at T/fj, — > 0, the integrals themselves are 



convergent in that limit, as they should. 
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